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ABSTRACT 

We derive the luminosity functions of elliptical galaxies, galaxy bulges, galaxy pseudo-bulges and 
galaxy discs from our structural catalogue of 10,095 galaxies. In addition we compute their associated 
luminosity and stellar mass densities. We show that spheroidal systems (elliptical galaxies and the 
bulges of disc galaxies) exhibit a strong color bimodality indicating two distinct types of spheroid which 
are separated by a core color of (u — r) ~ 2 mag. We argue that the similarity of the red elliptical 
and the red bulge luminosity functions supports our previous arguments that they share a common 
origin and surprisingly find that the same follows for the blue ellipticals and blue bulges, the latter of 
which we refer to as pseudo-bulges. In terms of the stellar mass budget we find that 58 ± 6 per cent is 
currently in the form of discs, 39 ± 6 per cent in the form of red spheroids (13 ± 4 per cent ellipticals, 
26±4 per cent bulges) and the remainder is in the form of blue spheroidal systems (~ 1.5 per cent blue 
ellipticals and ~ 1.5 per cent pseudo-bulges). In terms of galaxy formation we argue that our data on 
galaxy components strongly supports the notion of a two-stage formation process (spheroid first, disc 
later) but with the additional complexity of secular evolution occurring in quiescent discs giving rise 
to two distinct bulge types: genuine 'classical' bulges and pseudo-bulges. We therefore advocate that 
there are three significant structures underpinning galaxy evolution: classical spheroids (old); pseudo- 
bulges (young) and discs (intermediate) . The luminous nearby galaxy population is a mixture of these 
three structural types. The nature of the blue elliptical galaxies remains unclear but one possibility is 
that these constitute recently collapsed structures supporting the notion of mass-dependent spheroid 
formation with redshift. 

Subject headings: galaxies: spiral - galaxies: structure - galaxies: photometry - galaxies: fundamental 
parameters - ISM: dust, extinction 



1. INTRODUCTION 

In a recent paper (Driver et al. 2006) we demonstrated 
that galaxy bimodality is not just evident in color 
(see Strateva et al. 2001; Baldry et al. 2004) but also in 
the joint color-structure plane (see also Ball et al. 2006; 
Park et al. 2006; Choi et al. 2006; Conselice 2006). In 
that work we used SDSS photometry and single Sersic 
(1963; Graham & Driver 2005) profile fits to investigate 
the distribution of 10,095 relatively nearby luminous 
{Mb < —17 mag) galaxies in the color-Sersic index plane. 
The red peak is comprised of highly concentrated, high- 
Sersic index systems while the blue peak contains more 
diffuse, low-Sersic index systems. This is important as 
any movement from the blue peak to the red peak will 
require modifying the orbits, angular momentum and en- 
ergy of the entire stellar population. A simple inert pro- 
cess (e.g., exhaustion of the gas supply, stripping, etc.) 
could not achieve this, although a violent major merger 
event could (e.g., Barnes & Hernquist 1992). Perhaps 
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more importantly, when the population was segregated 
by Hubble type we found that while the early-type galax- 
ies (E/SOs, i.e., bulge dominated) lay almost exclusively 
in the red-concentrated peak, and the late-type spirals 
(Sd/Irr, i.e., disc dominated) in the blue-diffuse peak, 
the mid-type spirals (Sabc, i.e., bulge plus disc systems) 
straddled both peaks with no obvious sign of bimodality. 
We inferred from this that galaxy bimodality arises be- 
cause of the two component nature of galaxies and that 
spheroidal structures (i.e., ellipticals and bulges) will lie 
exclusively in the red-compact peak and discs in the blue- 
diffuse peak. As classical bulges lie within thin rotating 
disc systems this argues for early spheroid formation (via 
rapid merging or collapse) followed by a more quiescent 
phase in which the extended disc can form. To explore 
this further we have performed bulge-disc decomposition 
of all 10,095 galaxies in the Millennium Galaxy Cata- 
logue (see Allen et al. 2006) by fitting two-component 
Sersic-bulge plus exponential-disc models using GIM2D 
(Simard et al. 2002). In this Letter we report the lumi- 
nosity functions derived for various component samples 
(e.g. ellipticals, bulges, discs) and tabulate the associated 
luminosity and stellar mass densities for each component 
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class. 

Throughout this paper we assume a ACDM cosmology 
with Q m = 0.3, Q A = 0.7, and adopt h = H o /(100 km 
s -1 Mpc -1 ) for ease of comparison with other results. 

2. MGC COMPONENT LUMINOSITY FUNCTIONS 

The Millennium Galaxy Catalogue (MGC) is a deep 
(jUum = 26 mag arcsec~ 2 ), wide area (37.5 deg 2 ), £?-band 
imaging and redshift survey covering a 0.6 deg wide strip 
along the equatorial sky from lOh to 14h 50'. The MGC 
contains 10,095 galaxies down to B Mac = 20 mag, of which 
9,696 have redshift information. Full details of the imag- 
ing survey can be found in Liske et al. (2003), with the 
spectroscopic follow-up described by Driver et al. (2005; 
hereafter D05). In Allen et al. (2006) the bulge-disc de- 
compositions were reported for all 10,095 B MOC < 20 mag 
galaxies, and a final structural catalogue produced, con- 
sisting of bulge-only (Sersic profiles), disc-only (Sersic or 
exponential profiles) and bulge plus disc systems (Sersic 
plus exponential profiles). 

Luminosity functions are computed using the stan- 
dard step-wise maximum likelihood (SWML) estima- 
tor originally described by Efstathiou, Ellis & Peter- 
son (1988), with samples divided into bins of absolute 
magnitude. The MGC spectroscopic sample has a nom- 
inal Kron magnitude limit of B Mac = 20 mag. How- 
ever, here we will use the GIM2D total magnitudes 
(derived by integrating the light profiles to infinity) so 
that there is no longer a single limit that applies to 
the sample as a whole. To accommodate for this, each 
galaxy now has a unique magnitude limit, defined by: 
B lim = 20 + B MOC (Sersic) - B MCC (Kron). Following D05, 
we restrict our sample to galaxies in the redshift range 
0.013 < z < 0.18 and within carefully defined size and 
surface brightness boundaries (see D05 and Liske et al. 
2006 for full details). 

2.1. k+e corrections 

In D05 individual k-corrections were derived for each 
galaxy by comparing the total galaxy broad-band colors 
(uBgriz) to the 27 spectral templates given in Poggianti 
(1998) and identifying the best fitting spectrum. Having 
now separated the MGC galaxies into bulges and discs 
(see Allen et al. 2006) the global k-correction is not nec- 
essarily valid. However, approximately 50 per cent of 
our sample are best fit by one component profiles (i.e., 
bulge-only or disc-only galaxies) and for these systems 
we adopt the k-corrections as previously derived in D05. 
For the remaining, two-component systems we consider 
our component colors too coarse to be used to derive ro- 
bust individual k-corrections (as our decompositions are 
done in a single filter only). For the case of blue bulges 
and blue discs the global k-correction is likely to be ap- 
propriate for both (assumming that the blue bulge has 
formed from the disc). In the case of discs surrounding 
classical red bulges we note that for low-B/T systems 
the k-correction is likely to be appropriate for the discs 
but not the bulges. In the case of high-B/T systems 
we note that Peletier & Bacells (1996) report that such 
discs are typically redder. We therefore consider it ap- 
propriate to continue to adopt the global k-correction 
for both our discs and blue bulges. For the red bulges, 
however, we adopt the spectral template most frequently 
adopted by our single component red ellipticals (a Sa 15 



Gyr spectrum). This spectrum can be represented by 
a fourth order polynomial valid over the redshift range 
< z < 0.18 only:' 

k(z) = 3.86z + 12.13z 2 - 50.14z 3 . (1) 

We note that if we follow a similar procedure for the discs 
and blue bulges the implied characteristic turnover lumi- 
nosities (shown in Table 1) are systematically reduced 
by ~ 0.1 mag and the implied stellar masses reduced by 
~ 14 per cent. 

To model evolution we assume pure luminosity evolu- 
tion of the following form: 

L z=0 = L(l + z)-P, (2) 

where (3 is set to 0.75 for the global luminosity function 
(see Driver et al. 2005), 1 for blue components (discs, 
blue bulges, blue ellipticals), and 0.5 for red components 
(red bulges and red ellipticals). These values are based 
on the recent results reported in Zucca et al. (2006) 
for red and blue systems (their type 1 and 4 respec- 
tively). We do not model number evolution as the red- 
shift range is small and our merger rate estimates, based 
on dynamically close pairs within the MGC, are low (see 
De Propris et al. 2005 and De Propris et al. 2007). 

2.2. Luminosity functions 

Fig. 1 shows the luminosity distributions and Schechter 
function fits for our full galaxy sample (upper left), el- 
lipticals only (upper right; i.e. objects with B/T = 1 
after logical filtering, see Fig. 13 of Allen et al. 2006), 
discs (lower left) and bulges (lower right). Note that 
the Schechter function in the upper left differs from that 
shown in D05 because the magnitudes are now based 
on Sersic profiles integrated to infinity rather than Kron 
magnitudes. This difference is significant, resulting in a 
brighter M* value by about 0.1 mag but a comparable 
faint-end slope, a. Analysis of independent repeat ob- 
servations of ~ 700 galaxies suggests that our decompo- 
sitions are valid to good accuracy (AM BDLOE = ±0.1 mag 
and AM DISC = ±0.15 mag, see Allen et al. 2006) for com- 
ponentss with luminosities with Mb < —17 mag. Below 
this limit our decompositions become increasingly less re- 
liable and these data are shown with open symbols. The 
most striking result from Fig. 1 is the rapidly rising faint- 
end slope for the elliptical population. This was noted 
previously in Driver et al. (2006) and was shown to be 
due to contamination of the classical elliptical sample by 
low luminosity blue spheroids (see also Ellis et al. 2005). 

In Fig. 2 we show the color-structure plane defined by 
SDSS core (u — r) PSF color versus component Sersic 
index for the ellipticals (upper left) and galaxy bulges 
(upper right) . The bimodality of the ellipticals is strik- 
ing, with a blue and red population being apparent. The 
blue sample defines what we label blue ellipticals which 
were identified in Ellis et al. (2005) and quantified in 
Driver et al. (2006). A cut at (u — r) = 2 mag provides 
a clear division. The lower panels of Fig. 2 show the lu- 
minosity distributions and Schechter function fits for the 
blue and red samples. We see that red ellipticals and 
blue ellipticals follow markedly different trends and it is 
indeed the blue ellipticals which are responsible for the 
apparent upturn in the total elliptical galaxy luminosity 
function at very faint absolute magnitudes. We note that 
when the bulges are divided in the same manner, the blue 
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bulge luminosity function is very similar to that of the 
ellipticals, possibly indicating some common origin. 

It is tempting to associate the blue bulges with pseudo- 
bulges (see Kormendy & Kennicutt 2004) , which are be- 
lieved to arise from inner disc instabilities giving rise to 
a 'swelling' of the disc in the central region. As many of 
our blue bulges have Mb > — 17 mag, where our bulge- 
disc decompositions become unreliable, we cannot un- 
ambiguously confirm this population as pseudo-bulges 
but for the moment adopt this nomenclature for ease 
of discussion. The blue ellipticals remain somewhat in- 
triguing and appear to define a new class of object as 
previously noted by Ellis et al. (2005) and Driver et al. 
(2006; see also Ilbert et al. 2006 who identify them as a 
rapidly fading population). We are currently exploring 
these systems further (Ellis et al. 2007) and for the mo- 
ment simply flag them as interesting. From their distinct 
luminosity function it is clear they are predominantly low 
luminosity systems and could potentially represent the 
local counterparts to the luminous blue compact galax- 
ies studied by Guzman et al. (1997) and Phillips et al. 
(1997). 

Fig. 3 shows the final component luminosity functions 
with the red ellipticals and red bulges combined into a 
single red spheroid group and the blue ellipticals and blue 
bulges grouped together into a single blue spheroid class. 
The justification for this is the similarity in the shapes 
and ranges of the luminosity distributions from Fig. 2. 

2.3. Luminosity densities and stellar-mass densities 

Table 1 shows the Schechter function values for the lu- 
minosity function fits shown in Figs. 1, 2 & 3, along with 
their associated b j-band luminosity, jbj , and stellar mass 
densities, pm- These are derived using the following ex- 
pressions: 

j bj = 0*lO-°- 4(M «- Me) r(a + 2) (3) 

and 

N 

pm = J2 (0 4 /iv l )io (1 - 93(s - r)i - o - 79) io- o - 4(MB -- M e ) . 

i 

(4) 

The latter expression is first shown in Driver et al. (2006) 
and is based on the color to mass-to-light ratios given 
by Bell & de Jong (2001) which assume a Salpeter-'lite' 
IMF. Note that the (g — r) colour for each galaxy was 
obtained by matching the MGC to the Sloan Digital Sky 
Survey first data release (Abazajian et al. 2003). 

3. DISCUSSION 

From Table 1 we see that 58 ± 6 per cent of the stellar 
mass is in the form of galaxy discs, 13 ± 4 per cent in 
red elliptical galaxies and 26 ± 4 per cent in classical red 
bulges. Previously it has been reported (Bell et al. 2003; 
Baldry et al. 2004; Driver et al. 2006) that 54 - 60 per 



cent of the stellar mass lies in the red peak dominated by 
early-type galaxies. As we have now separated the early- 
type galaxies into bulges and discs one expects that the 
spheroid (elliptical+bulge) stellar mass density should 
be lower than the early-type stellar mass density, which 
of course includes the associated disc components of the 
lenticular galaxies. Examining the color-structure plane 
(Fig. 2 upper panels) for ellipticals and bulges we see 
that the red populations (of each type) lie in the same 
location. This is indicative of a shared origin for red 
ellipticals and classical red bulges. When combined the 
red spheroids account for 39 ± 6 per cent of the stellar 
mass density. Hence the bulk of the stellar mass (96 
per cent) exists in the two classical structures originally 
defined by de Vaucouleurs (1959). 

The remaining 3 per cent lie in the form of blue ellip- 
ticals (~ 1.5 per cent) and blue bulges (~ 1.5 per cent). 
These latter two populations can therefore be considered 
minor from a cosmological perspective. Furthermore, as 
the blue bulges are likely to represent either difficulties 
in the decomposition (e.g., bars) or pseudo-bulges (disc 
swelling), their stellar mass could arguably be added, 
in either case, to that of the galaxy discs (i.e., ~ 61 
per cent). The nature of the blue ellipticals (previously 
dubbed blue spheroids) remains uncertain but they ap- 
pear to constitute a very small fraction of the stellar 
mass budget, although we must note the near divergent 
faint-end slopes. The primary conclusion then is that 
the stellar mass is mainly divided between two distinct 
structures: blue 2D discs and red 3D spheroids. 

Finally we note that the results presented in this Letter 
are based on B-band data and therefore susceptible to 
dust attenuation (see Shao et al. 2006). Bell & de Jong 
(2001) argue that the effect on an individual galaxy's stel- 
lar mass is less than one might expect (see their fig. 1) 
because the observed decrease in total luminosity is off- 
set by the increased stellar mass-to-light ratio inferred 
from the redenned colors, therefore yielding comparable 
final stellar masses. We explore this in detail in Driver 
et al. (2007) and note that while the masses are mod- 
ified somewhat the final stellar mass breakdown is not 
dramatically altered. 



The Millennium Galaxy Catalogue consists of imag- 
ing data from the Isaac Newton Telescope and spectro- 
scopic data from the Anglo Australian Telescope, the 
ANU 2.3m, the ESO New Technology Telescope, the 
Telescopio Nazionale Galileo, and the Gemini Telescope. 
The survey has been supported through grants from 
the Particle Physics and Astronomy Research Coun- 
cil (UK) and the Australian Research Council (AUS). 
The data and data products are publicly available from 
http://www.eso.org/~jliskc/mgc/ or on request from 
J. Liske or S.P. Driver. 
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TABLE 1 



bamplc 


M B — o log n 
(mag) 


a 


•P 

(10 A Mpc 5 (0.5 mag) 1 ) 


jb,, 

(10 8 hL Q Mpc -3 ) 
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AT 
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Fig. 1. — B-band luminosity functions derived using our variable limit SWML method for the global sample (top left), discs (bottom 
left), ellipticals (top right), and bulges (bottom right). In each case the derived data (squares), and fitted Schcchtcr functions (solid lines) 
are shown. The Schechter fit for the global sample uses all data points but for the component luminosity functions (discs, ellipticals, and 
bulges) we only use data points with Mg < —17 mag (see Section 2.3). 
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Fig. 2. — The top row shows the (u — r) — n distribution (dots) for ellipticals (left), and bulges (right) with Mb < — 17 mag. The contours 
show the volume-corrected luminosity density for these objects. The natural division between the red and blue populations appears to 
occur at (u — r) = 2 mag. The bottom row shows the B-band luminosity functions for ellipticals (left), and bulges (right), split into red, 
(u — r) > 2 mag, components where data points are shown by squares, and blue, (u — r) < 2 mag, components where data points are shown 
by triangles. 
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Fig. 3. — Derived data and Schechter fits for B-band luminosity functions for our three final structural components: discs (circles), 
red bulges and ellipticals (squares), and blue (pseudo-)bulges and blue ellipticals (triangles). The Schechter fits only use data points with 
Mb < — 17 mag. The thicker black line shows the global B-band luminosity function. 



